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The rich phenomenology of collective neutrino oscillations has been studied only in one-dimen- 
sional or spherically symmetric systems. Motivated by the non-spherical example of coalescing 
neutron stars, presumably the central engines of short gamma-ray bursts, we use the Liouville 
equation to formulate the problem for general source geometries. Assuming the neutrino ensemble 
displays self-maintained coherence, the problem once more becomes effectively one-dimensional along 
the streamlines of the overall neutrino flux. This approach for the first time provides a formal 
definition of the "single-angle approximation" frequently used for supernova neutrinos and allows 
for a natural generalization to non-spherical geometries. We study the explicit example of a disk- 
shaped source as a proxy for coalescing neutron stars. 
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I. INTRODUCTION 

Flavor transformations caused by neutrino mixing de- 
pend on the matter background and on the neutrino 
fluxes themselves: neutrino-neutrino interactions pro- 
vide a nonlinear term in the equations of motion [l|, Q 
that gives rise to collective flavor transformation phenom- 
ena 0, i, i, i, 0, i, i, M, El, [H, m • The neutrino den- 
sity needs to be so large that a typical neutrino-neutrino 
interaction energy fi is comparable to the vacuum oscil- 
lation frequency u ~ Am^/2E. Only recently has it 
been fully appreciated that this condition is sufficient 
even if a dense background of ordinary matter provides a 
much larger interaction energy so that naively neutrino- 
neutrino interactions would seem negligible [ij, • Fol- 
lowing this crucial insight, nonlinear oscillation phenom- 
ena in the supernova (SN) context have been studied over 
the past t wo y ears in a long series of p apers fisl. [igI . [itI . 

m, m, HO, mills, m, m, m, m, 113, m M mm- 

One striking effect is "self-maintained coherence" d, 
0, H, @ ■ Different neutrino modes have different vacuum 
oscillation frequencies uj = Am^/2E, but with strong 
neutrino-neutrino interactions they "stick together" and 
oscillate as a single mode characterized by the "synchro- 
nized oscillation frequency" Wsync = {Arn^ /2E). This 
can lead to all modes going through an MSW resonance 
together, the "collective MSW-like transition" or "syn- 
chronized MSW effect" d, i, US, HI, HI] . 

More interesting still are collective phenomena driven 
by the decrease of the neutrino flux with distance from 
the source. The adiabatic transition from a dense to a 
dilute neutrino gas produces step-like spectral features 
where the spectrum sharply splits into parts of different 
flavor transformation, so-called "step- wise spectral swap- 
ping" or "spectral splits" [H, M, 0, 113, Hj, 113, IH, Ha, 
[29j . Spectral splits can result from a preceding collective 
MSW effect ("MSW prepared spectral split") or from 
neutrino-neutrino interactions alone. 

The latter case depends on an unusual form of non- 



equilibrium among neutrino flavors where one has an ex- 
cess of flavor pairs, say i^ei^e, over the other flavors. For 
neutrinos streaming off a SN core one indeed expects a 
hierarchy of number fluxes F,^^ > Fp^ > F^^^^^ = Fp^^p^. 
Therefore, one can have collective transformations of the 
form Ve^e Vx^x, where x stands for some suitable com- 
bination of ^ and r neutrinos. These "collective pair 
transformations" do not violate any conservation law and 
thus can be catalyzed even by a very small mixing an- 
gle. Fp^ can completely swap with Fp^ whereas the larger 
Fp^ converts only to the extent allowed by flavor-lepton 
conservation, but in a step-like spectral form. Both the 
complete conversion of and the split in the u,, spec- 
trum provide signatures for the inverted neutrino hierar- 
chy even for an extremely small 13-mixing angle (23l.l30j. 

For non-isotropic enviroments, multi-angle effects may 
play an important role. The term "multi-angle effects" 
actually refers to two different issues. One is that the 
weak interaction potential between two relativistic parti- 
cles is proportional to (1 — cos^), where is their rel- 
ative angle of propagation. One usually considers an 
isotropic background of ordinary matter so that cos 9 av- 
erages to zero. The same is true in an isotropic neu- 
trino gas for the neutrino-neutrino term. The second 
issue is the "multi-angle instability" . Neutrinos arriving 
from different points on the source belong to different 
angular modes , wh ich may decohere kinematically in fla- 
vor space ini, [13, El, El Hii- This effect can be self- 
induced in the sense that a very small initial anisotropy 
is enough to trigger an exponential runaway, for exam- 
ple in a gas consisting of equal densities of neutrinos 
and antineutrinos (l7j . Systems consisting of very few 
angular modes can show a two-stream or multi-stream 
instabihty [ll|, E3, El- 

On the other hand, numerical 
studies show that systems consisting of many angular 
modes and with a sufflcient neutrino-antineutrino asym- 
metry do not show a multi-angle instability but rather 
show self-maintained coherence among different angular 
modes [l^ HH, |13|- The SN neutrino flux parameters 
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seem to be such that the multi-angle instability plays no 
role in practice. Based on this assumption, most of the 
SN studies have used the "single-angle approximation" , 
where all angular modes are assumed to have the same 
behavior. 

Multi-angle effects related to the (1 — cos 6*) structure 
of the neutrino-neutrino term are unavoidable for an ex- 
tended source radiating neutrinos into space because the 
emitted neutrino flux cannot form an isotropic gas. How- 
ever, these effects also occur, and are easier to study the- 
oretically, in a homogeneous system evolving in time that 
has a non-isotropic angular distribution of neutrinos. 

In practice, however, one usually deals with station- 
ary systems where one asks for the spatial variation of 
a neutrino ensemble as a function of distance from the 
source. Even if the neutrino-neutrino interaction were 
isotropic, we still would have geometric multi-angle ef- 
fects because neutrinos reaching a certain point from an 
extended source have traveled on different trajectories. 
Even the simple case of an infinite radiating plane is not 
trivial. Here the direction perpendicular to the plane is 
the only direction in which the overall neutrino ensem- 
ble can show any spatial variation. Even if all neutrinos 
have the same energy and thus oscillate with the same 
frequency along their trajectories, the projection on the 
direction perpendicular to the plane yields different effec- 
tive oscillation frequencies and thus kincmatical decohcr- 
ence. Neutrino-neutrino effects can synchronize different 
angular modes so that a sufficiently dense neutrino gas 
will not show this form of multi-angle decoherence. On 
the contrary, all angular modes will vary with the same 
oscillation length as a function of distance from the plane. 
A similar description applies to a spherical source where 
one asks for the variation of all angular modes along the 
radial direction. 

The single-angle treatment of SN neutrino oscillations 
amounts to the assumption of self-maintained coherence 
among angular modes, although it has never been explic- 
itly expressed in this form. Some authors assumed that 
all angular modes oscillate as the radial one How- 
ever, in this case the neutrino-neutrino interaction van- 
ishes because of the (1 — cos 6') factor, so it was necessary 
to assume a certain average of the neutrino-neutrino in- 
teraction strength. Other authors represented all angular 
modes by a single angular mode radiated at 45° relative 
to the radial direction and then used a neutrino-neutrino 
interaction strength consistent with this assumption [27j . 
This implementation of the single-angle approximation 
has the advantage that one can use the same numeri- 
cal code as for multi- angle simulations, simply restrict- 
ing oneself to a single angular bin. In yet other cases 
the system was modeled as a homogeneous and isotropic 
gas that evolves in time, assuming a time variation of the 
neutrino density that mimics the radial variation in the 
spherical case. 

One of our goals is to show that the single-angle treat- 
ment can be formulated self-consistently. The assump- 
tion that all angular modes evolve the same in flavor 



space provides a unique concept of what is meant by 
"single-angle behavior." This is straightforward in the 
systems described so far where symmetry dictates that 
the spatial variation is only along a certain direction, ef- 
fectively reducing the problem to one dimension. We are 
really motivated, however, by more general geometries 
where no special direction is singled out by symmetry. 
In particular, we are interested in the case of coalescing 
neutron stars that may form the inner engines of short 
gamma-ray bursts [s^ . 

The accretion torus or disk formed during neutron star 
coalescence is a neutrino source comparable to a SN core. 
However, the torus is less dense and not efficient at pro- 
ducing and Vr- Therefore, the torus is a source for 
a dominant Ve^e pair flux which is thought to produce 
an e^e~ pair plasma, thus powering short gamma-ray 
bursts. The annihilation cross section for v^De — > e^e~ 
is much larger than that for Vx^x e~^e~, so the neu- 
trino flavor composition strongly influences the number 
of e+e^ pairs produced. Therefore, one may ask if col- 
lective pair conversions occur in this environment close 
enough to the source to modify the energy transfer to 
the e"^e~ plasma, and hence affect the strength of the 
gamma ray burst. 

For coalescing neutron stars one expects a flux hierar- 
chy Fp^ > F^^ > F^^ = Fp^ , which differs from the SN 
case because the matter leptonizes when neutrons con- 
vert to protons, in contrast to the deleptonization of a 
SN core. The asymmetry between F,^^ and Fp^ could be 
enough to prevent multi-angle decoherence so that sim- 
ilar collective effects as in the SN environment are con- 
ceivable. However, even granting this assumption, it is 
not straightforward how to implement something like a 
single-angle approximation in this context because it is 
not obvious how one should picture self-maintained co- 
herence. 

The purpose of our paper is to formulate the mean- 
ing of self-maintained coherence for general source ge- 
ometries and study its implications. We flnd that the 
flavor variation reduces to a quasi one-dimensional prob- 
lem along the streamlines of the total neutrino flux. The 
main difference between the general case and the radi- 
ating plane or sphere is that the streamlines are typi- 
cally curved, at least close to the source, so that self- 
maintained coherence applies to flavor oscillations along 
these curved streamlines. 

We begin in Sec. |TT] with the general equations of mo- 
tion for the neutrino matrices in flavor space. In Sec. Ill Bl 
we formulate the collective equations for neutrinos only 
(no antineutrinos) and consider only synchronized oscil- 
lations. Sec. Ill CI shows the existence of streamlines and 
gives a prescription for calculating the flavor evolution 
along them. In Sec. IIIII we solve the problem for several 
geometries. In Sec. IIVI we study the generalization to a 
mixed system of neutrinos and antineutrinos where col- 
lective pair transformations are possible. In Sec. |V] we 
consider an explicit example for a disk source with pa- 
rameters inspired by numerical simulations of coalescing 
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neutron stars. Wc summarize our conclusions in Sec. lVIl 

II. FORMALISM FOR AN ARBITRARY 
SOURCE GEOMETRY 

A. General framework 

A homogeneous ensemble of mixed neutrinos can be 
described by "matrices of density" pp, which really are 
matrices of occupation numbers, for each momentum 
mode p 0, m, m, HI]. If and a^^p are the cre- 
ation and annihilation operators of a neutrino in the mass 
eigenstate i of momentum p, we have (pp)y cx; (a|ai)p 
so that the diagonal entries of Qp are the usual occupa- 
tion numbers (expectation values of number operators), 
whereas the off-diagonal elements encode the phase re- 
lations that allow one to follow flavor oscillations. Such 
a description assumes that higher-order correlations be- 
yond field bilinears play no role, probably a good approx- 
imation for neutrinos produced from essentially thermal 
sources such as the early-universe plasma or a SN core. 

Antineutrinos arc described in an analogous way by 
{Qp)ij = {0'iO,j)p- Note that we always use overbars 
to characterize antiparticle quantities. The order of fla- 
vor indices was deliberately interchanged on the r.h.s. 
so that the matrices Pp and Qp transform in the same 
way in flavor space [4]. In this way one can, for ex- 
ample, write the overall neutrino current in the simple 
form / dp-(;p(pp — Qp) where here and henceforth dp 
stands for d'^p/(27r)^. While some authors prefer the 
seemingly more intuitive equal order of flavor indices for 
neutrinos and antineutrinos, in that convention the cur- 
rent would be / dp?;p(ep — p*) and generally the equa- 
tions will involve both the matrices of density and their 
complex conjugates. (In a truly field-theoretic derivation 
there is no ambiguity about the relative structure of neu- 
trino and antineutrino matrices. An analogous example 
for matrices of occupation numbers is provided by the 
kinetic treatment of the quark-gluon plasma where the 
quark distribution functions are 3x3 matrices in color 
space [13. The matrices for quarks and the ones for 
antiquarks transform equally under a color gauge trans- 
formation.) 

The matrices Qp and Qp depend on time. Their evolu- 
tion is governed by the Boltzmann collision equation 

Step = -iPp,(?p] +£'p|^o„ (1) 

and a similar equation for ^p. The collision term will be 
of no further concern in our paper because wc only study 
freely streaming neutrinos. Further, [•, •] is a commutator 
and f2p is the matrix of oscillation frequencies. In vacuum 
we have f2p = M^/2|p| with M the neutrino mass matrix. 
In general Qp also depends on the background medium 
and notably on the presence of other neutrinos. Here 
and henceforth the ultrarelativistic limit for neutrinos is 
assumed. In other words, it is assumed that the difference 
between the neutrino energy Ep and momentum p ~ 



IpI is irrelevant except for oscillations: the only relevant 
difference between energy and momentum is captured by 
the matrix of oscillation frequencies. 

Up to this point we have considered a homogeneous 
ensemble evolving in time, largely relevant for neutrino 
oscillations in the early universe for which this formal- 
ism was originally developed. However, for a realistic 
representation of radiating objects such as supernovae 
or coalescing neutron stars we need to include spatial 
transport phenomena. To this end one introduces space- 
dependent occupation numbers or Wigncr functions Qp.x- 
The quantum-mechanical uncertainty between location 
and momentum implies that this construction is only 
valid in the limit where spatial variations of the ensemble 
are weak on the microscopic length scales defined by the 
particles' typical Compton wavelengths. The left hand 
side of the Boltzmann collision equation now turns into 
the usual Liouville operator [3^ 

Stgp.x + Vp • Vx t'p^x + P • Vp Pp.x- (2) 

The transition to the Liouville operator may seem obvi- 
ous, but making it conceptually precise requires a long 
argument (STj . The first term represents an explicit time 
dependence, the second a drift term caused by the parti- 
cles' free streaming, and the third the effect of external 
macroscopic forces, for example gravitational deflection. 
Henceforth we shall neglect external forces and only re- 
tain the drift term. 

In our ultrarelativistic approximation the modulus of 
Vp is the speed of light. Of course, for propagation over 
very large distances this may be a bad approximation 
when time-of-flight effects play a role. In this case the 
drift term would have a more complicated structure be- 
cause the velocity is then also a nontrivial matrix in flavor 
space 0. 

We shall focus on stationary problems. Even neutrino 
emission from a SN or coalescing neutron stars fall in this 
category unless there are fast variations of the source. 
In other words, we shall ignore a possible explicit time 
dependence of gp x so that flnally we arrive at 

Vp • Vx Qp.x = -i[f2p,x, gp,x] , (3) 

where the matrix of oscillation frequencies in general also 
depends on space because of the influence of matter and 
other neutrinos. 

Equation ^ could have been guessed starting from the 
usual treatment of neutrino oscillations. The neutrino 
density matrix describing a stationary beam in the z- 
direction follows dzQE,z = —'AP'E.z,Qe,z\- Moreover, in 
our ultrarelativistic approximation E = p and |vp| = 1. 
Equation ([5]) then simply amounts to a collection of many 
beams in different directions originating from different 
source points. 

We finally spell out Eq. ^ explicitly under the as- 
sumption that the only nontrivial ingredients consist of 
neutrino-neutrino refractive effects. 
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iVp-Vxffp.x = +^ [M^, gp^x] + \/2Gf J dq(l - Vq • Vp)[(£(q,x - Pq.x), Pp^x] , 
iVp-Vx^p,x = [M^, ^p.x] + %/2Gf J dq(l - Vq • Vp)[(£lq,x - eq,x), ^p,x] • 



(4) 



r 



The only difference between the neutrino and antineu- 
trino equations is the sign of the vacuum osciUation term. 
This structure is a consequence of defining the antineu- 
trino matrices with reversed flavor indices. 



B. EOMs with self-maintained coherence 

The nonlinear equations of motion (EOMs) Eq. (|4]) 
simplify considerably if self-maintained coherence occurs 
in a dense neutrino gas and all modes can be assumed to 
evolve in the same way. For the moment we restrict our- 
selves to a source radiating only neutrinos (no antineu- 
trinos). What we mean by self- maintained coherence in 
this case is that at a given location all neutrino modes 
are aligned with each other, assuming they were aligned 
at the source. That is, 



^?p,x Px /p.x 



(5) 



Here, /p.x ~ Tr(gip.x) is a scalar occupation number, 
summed over all flavors, while for N flavors Px is a iV x 
matrix normalized as Tr(Px) = 1. We also define the 
local neutrino density and flux as 



where 



"x = y dp/p,x, 

Fx = y"dpvp/p,x = (v)x'^3< 



^^^^^/dpQp,x/p,: 



/dp/p,: 



(6) 



(7) 



is the momentum average of a quantity Q at location x 
with respect to the distribution function /p.x- For future 
convenience we also define 



(8) 



a unit vector along the direction of flux Fx, or equiva- 
lently, along the average velocity (v)x. 

In this language, the EOM Eq. Q becomes 

Vp-Vx/p,xPx = -^[M2,Px]/p,x, (9) 

where the nonlinear term has disappeared because it in- 
volves the commutator [Px, Px] = 0. The only impact of 



the neutrino-neutrino term is to glue the different modes 
together, but afterwards it no longer appears in the equa- 
tion. 

The evolution now factorizes into the free streaming of 
the overall neutrino flux and the oscillation of the com- 
mon flavor matrix Px- The trace of Eq. ^ provides 



Vx/p,: 



0, 



(10) 



allowing us to determine the overall neutrino density 
and flux Fx at any point, assuming it is given at the 



source surface. With Vp • Vx /p,x = 
grating Eq. pUj) over all modes yields 

Vx • Fx = , 



p ./p,X: 



mte- 



(11) 



expressing the absence of neutrino sources or sinks. 

The EOM for Px is found by integrating Eq. (O over 
all modes, 



Fx P. 



2 L 



(12) 



The l.h.s. expands as (Vx • Fx) Px + Fx • (VxPx) where 
the first term disappears because of Eq. pl|) so that 



VxP. 



dp 



/p 



Dividing both sides by |Fx|; we get 

Fx • VxPx = -i [f2x, Px] 



where 



^ _ Jdpp-\fp, 
2 |/dpvp/p,x 
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l(v) 



(13) 



(14) 



(15) 



Equation p4)) resembles the EOM in the single-angle ap- 
proximation, with fix as the synchronized matrix of os- 
cillation frequencies. Note that {p~^} is independent of 
location if the energy spectrum of neutrinos is assumed 
to be the same everywhere. 



C. Streamlines 

Equation is a partial differential equation (PDF) 
for the matrix Px. It can be reduced to a set of ordinary 
differential equations (ODEs) 



dx 
dJ 
dPx 

ds 



-i [f^x, Px] , 



(16) 
(17) 
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where s is a parameter along the "characteristic Une," 
or "streamhne," defined by Eq. ([16]). Since Fx is unique 
at each location, the streamlines do not intersect each 
other. Along each streamline, the differential equation 
[Eq. pT)) ] for the matrix Px is a set of coupled, but linear 
ODEs and hence can be solved easily and uniquely, given 
the boundary conditions. 

As an illustration, we now specialize to a two-flavor 
system, where the coupled ODEs in Eq. p?]) reduce to a 
single ODE and may be interpreted as the evolution of 
a single phase. As uaual, we express a Hcrmitian 2x2 
matrix A in terms of a polarization vector A by virtue of 
A = i (Tr A + A - a) with a the vector of Pauli matrices.^ 
The equation for synchronized oscillations then turns into 
the usual precession formula 



VxPx = CJx S X Px 



(18) 



where Px is the polarization vector corresponding to 
and uj^B the one corresponding to Qx- Here P is a 
unit vector by definition, which is determined by and 
hence is independent of location. The effective synchro- 
nized oscillation frequency is 



1^ 



(19) 



where wq is the synchronized oscillation frequency far 
away from the source, where the details of the source 
geometry do not matter. All the effects of the source 
geometry are thus captured by the quantity |(v)x|. 

The problem becomes even simpler if we observe that 
the magnitude of Px and its projection on B are con- 
served. We can then write Px in terms of its components 
along and perpendicular to B as 

Px = {P± cos ax, P± sin a^, Pj, ) . (20) 
Equation then reduces to a scalar equation 



Vx flx 



(21) 



In other words, Px is fully specified by a real scalar num- 
ber flx that tells us the phase of the revolution of the 
polarization vector around B. Equations pB]) and pT]) 
in this case simplify to 



dx 
d^ 
dflj 
"dl 



l(v)> 



(22) 
(23) 



As before, the streamlines are given by Eq. ([22]) and 
the evolution is obtained by solving the single ODE in 



Eq. p3|) . This has been possible because of a simple 
parametrization of Px in Eq. (j20p in two flavors. 

In the case of three flavors, one can express a Hermitian 
3x3 matrix in terms of the Gell-Mann matrices [1^. 
Equation psp then retains the same form, except that 
the "cross product" is now understood to be {A x B)i = 
fijkAjBk with fijk the structure constants for SU(3). 
The ODEs in Eq. p7|) are coupled, but linear, and can 
be solved along the streamlines. 

The general strategy for calculating the neutrino flavor 
evolution is therefore the following, (i) Find the stream- 
lines for the problem by solving Eq. p6[) with the given 
neutrino source, (m) Solve for the phase (or for the 
matrix Px) by solving the ODEs [Eq. p?|) ] along each 
streamline. This leads to the determination of ^ ev- 
erywhere. 

An important consequence of the existence of stream- 
lines is that the evolution along any streamline is in- 
dependent of the other streamlines. Moreover, each 
streamline is determined solely by the flux Fx along it. 
Thus the neutrino flavor content at any point in space 
is determined simply by identifying the streamline that 
passes through that point and evolving the phase 
along the streamline. This effectively reduces any multi- 
dimensional problem to one dimension. 



III. FLAVOR EVOLUTIONS FOR VARIOUS 
SOURCE GEOMETRIES 

We now calculate the flavor evolution of neutrinos in 
several physical situations, ranging from the trivial to the 
complicated. The crucial quantity to calculate in each 
scenario is (v)x, which captures the effect of the source 
geometry through Eq. p9|) . 



A. Infinite Plane 

For a homogeneous plane that radiates neutrinos iso- 
tropically, the streamlines are straight lines normal to 
the surface. Therefore, the oscillation phase a depends 
only on the distance z from the plane and the surfaces 
of equal phase are planes parallel to the radiating sur- 
face. The average velocity of all neutrinos away from 
the plane is half the speed of light, (v)x = every- 
where. Therefore, Wx = 2a;o represents the flavor oscil- 
lations along the direction perpendicular to the plane. 
This gives a = 2ujo z + ao, where aq may be chosen to 
be an arbitrary real number. This, along with the initial 
conditions, completely determines the flavor state of the 
neutrinos at any point in space. 



B. Sphere 

A sans-serif letter such as A indicates a matrix in flavor space, a 

letter with an arrow such as A indicates a vector in flavor space. For a spherical source. Symmetry dictates that the 

and a bold-faced letter such as A a vector in coordinate space. Streamlines are Straight lines following the radial distance 
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FIG. 1: Synchronization frequency ujr for an emitting sphere 
(radius R) as a function of distance r. 



from the center and that the surfaces of equal phase are 
spherical shells around the source. Therefore spherical 
symmetry dictates that the average velocity (v)x = {vr)'r 
is along the radius. All that remains to be determined 
is (vr)-, the average radial velocity of all neutrinos at dis- 
tance r. 

To this end we note that all modes that are radiated 
at the neutrino sphere with the same zenith angle dfi 
relative to the radial direction behave identically. Their 
radial velocity at the neutrino sphere is vji = cos-dn, 
so it is useful to classify all modes by their vji. Simple 
geometry tells us that R sin'dii ~ rsiniJ^, where i?r is 
the zenith angle of a given neutrino trajectory relative to 
the radial direction at a distance r. Therefore, the radial 
velocity at distance r of a mode that has radial velocity 
VR at the neutrino sphere is 



cos -dr 



1-^(1-4) 



(24) 



Given the all- flavor phase-space density f{vR,R) of a 
mode VR at the neutrino sphere, we can determine its 
value f{vR, r) at another distance from flux conservation. 
In spherical coordinates we have dr[r'^Vrf{vR,r)] = so 
that f{vR,r) cx {R/rfivR/vr)- 

If we assume isotropic emission at the neutrino sphere, 
the modes are uniformly distributed in the variable < 
cos < 1 or < Wfl < 1 , so that 



S^dvR Vr f{vR,r) _ Jq dvR vr 
Jo <ivRf{vR,r) J^dyRVR/vr 



(25) 



Inserting the expression Eq. ()24p for Vr and performing 
the integral we find 




(26) 



iden- 



At the neutrino sphere (r = R) one finds (vr) = ^, 
tical to the radiating plane, whereas at large distances 
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FIG. 2; Streamlines for the disk source, combined with con- 
tours for various quantities: (a) Neutrino flux |Fx|. (b) Aver- 
age velocity |(v)x|. (c) Effective neutrino- neutrino interaction 
strength /ix. 



r 3> i? all modes move essentially in the radial direction 
and (vr) 1. Therefore, near the neutrino sphere, the 
neutrino flavor evolves along the radial direction with fre- 
quency uJr = 2uJo, whereas at large distances it evolves 
with This case is represented in Fig. [TJ The 

decrease of uJr to its asymptotic value is surprisingly fast. 
As usual, the phase is given by a = J uj^dr + ao, and this 
along with initial conditions specifics the flavour content 
of neutrinos all over space. 
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FIG. 3: Behavior as a function of z for = (left panels) and as a function of O at infinite distance (right panels) of the 
following quantities for the disk source: (a) Neutrino flux |Fx|. (b) Average velocity |(v)xl. (c) Effective neutrino-neutrino 
interaction strength /ix- 



C. Disk 

Perhaps the simplest example where the streamlines 
are curved is the case of a radiating circular disk. This 
model approximately mimics the accretion disk formed 
by two coalescing neutron stars. The disk radius is taken 
to be i?, setting the geometric scale for the system. We 
assume that the disk radiates homogeneously a flux _Fo 
throughout its surface and that at each point the emission 
is identical to a blackbody source. 

Very close to the disk the physical situation is identical 
to an infinite radiating plane. The streamlines are per- 
pendicular to the disk with a flux Fq. Viewed from a large 
distance, on the other hand, the flux must be directed 
away from what looks like a point-like source, so that the 
streamlines are radially outward from the source. Since 
the source is assumed to behave like a "blackbody" , the 
neutrino flux coming from the source varies with cos 
where Q is the zenith angle relative to the disk axis. This 
is simply because we are seeing a blackbody surface at an 
angle so that the total flux is reduced by the projected 
size of the source. On the other hand, the total flux pass- 
ing through a spherical surface at distance r is conserved. 



Therefore, at a distance r ^ R wc find 

|Fx| =Fo(i?/r)2 cose. (27) 

Therefore, wc know the flux and its direction directly at 
the source and at a large distance. 

We calculate the flux Fx numerically and find that the 
streamlines are well represented by the hyperbolae 



sin^ 9 cos^ 9 



(28) 



Here we use coordinates x along the disk radius and z 
along the disk axis. We have also calculated the mag- 
nitude of fiux I Fx I, the magnitude of average velocity 
|(v)x| = |Fx|/?^xj and the effective neutrino- neutrino in- 
teraction strength fi^ defined later in Eq. (|3T|l . Wc show 
contours of these quantities in Fig.[21 Moreover, in Fig. [3] 
we show the behaviour of these quantities as a function 
of z for 9 = and their asymptotic variation as a func- 
tion of 9 at a large distance from the disk. Details of the 
numerical calculations are given in Appendix |^ 
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IV. NEUTRINOS AND ANTINEUTRINOS 
A. Single Mode Equations 

We now extend the previous considerations to the more 
interesting situation where both neutrinos and antineu- 
trinos stream from the source. The main difference is 
that we now assume that neutrinos and antincutrinos 



each form two separate classes of modes that oscillate 
together, i.e., we assume that the behavior of all neu- 
trino modes is well represented by a single matrix Px 
and that of the antineutrino modes by a matrix Px such 
that £>p,x = Px/p,x and pp^x = Px.fp,x. 

Integrating the EOMs of Eq. ([4]) over all modes we 
arrive at 



iFx • Vx Px = +[^0, Px]nx + V2Gf ([(rixPx - nxPx),nxPx] - [(FxPx - FxPx),FxPx]) 

iFx- VxPx = -[fio,Px]r^x + %/2GF([(nxPx-nxPx),?lxPx] - [(FxPx-FxPx),FxPx]) 

I 



(29) 



If it were not for the nonlinear term, of course, we would 
have synchronized oscillations separately for neutrinos 
and antincutrinos where f2o is the oscillation matrix av- 
eraged over the neutrino spectrum whereas Q.q is the one 
averaged over the antineutrino spectrum. 

In order to consider the simplest nontrivial case we as- 
sume further that the energy spectra for neutrinos and 
antincutrinos are the same as well as their angular emis- 
sion characteristics. This implies fip = ^o- Moreover, 
we assume that the total flux of antincutrinos is a < 1 
times the flux of neutrinos. It simplifies the equations if 
we absorb this global factor in the normalization of the 
matrix Px so that Tr(Px) = a whereas Tr[Px) = 1. With 
these assumptions we have = and Fx = Fx. With 
(v)x = Fx/^x = (v)x we then find 

iFx-VxPx = +[fix,Px] +Aix[(Px-Px),Px], 
iFx-VxPx = -[fix,Px] +A^x[(Px-Px),Px].(30) 

The oscillation matrix is f2x = Qo/|(v)x| whereas /ix is 
the effective neutrino-neutrino interaction strength to be 
discussed below. 

If neutrinos and antincutrinos have the same stream- 
lines, we find a unique extension of the idea of self- 
maintained coherence to a source with general geome- 
try. Therefore, one may expect that the solutions along 
given streamlines depend on the source geometry primar- 
ily through the variation of /i along the streamlines. 

In the SN context, the evolution of the neutrino en- 
semble driven by the decrease of /i tends to be essentially 
adiabatic so that the exact variation of fi along a stream- 
line should not matter much for the overall picture. This 
suggests that we may be allowed to borrow the insights 
gained in the SN context directly to more general cases. 

Even if neutrinos and antincutrinos do not follow ex- 
actly the same streamlines, if the lateral variation of the 
solution is smooth, one may speculate that the system 
finds some suitable average and that the final solution 
could still be similar to the simple case. 



B. Effective Neutrino-Neutrino Interaction 
Strength 

In the limit of self-maintained coherence discussed in 
the previous section, we are led to define an effective 
neutrino-neutrino interaction strength 



Aix = V2Gf|Fx 
= x/2Gf|Fx 



<-|Fx 
iFxP 
1 



- 1 



(31) 



In principle, the same coefficient would have appeared 
in the neutrino-only case, but there is was accompanied 
by a vanishing commutator so that it did not appear ex- 
plicitly. The definition of /ix is unique except for the 
normalization that depends on our choice Tr(Px) = 1 
and Tr(Px) = a. In our normalization, /Xx has the in- 
terpretation of the average neutrino-neutrino interaction 
energy and corresponds to the quantity v^Gpn-e in the 
ordinary matter effect. 

Equation (j3ip implies some general features of the vari- 
ation of /i at large distance r from an extended but finite 
source. (We always consider sources that are not point 
like or else there would be no neutrino-neutrino effects, 
but that are finite, in contrast to the infinite plane men- 
tioned earlier, or else there would be no systematic de- 
crease of fi with distance.) The flux factor |Fx| provides 
a trivial scaling from the geometric fiux dilution, as 
shown in Fig. [31 The "coUinearity factor" 



Gx 



l(v) 



1 



(32) 



arises from the 1— cos 6 structure of the neutrino-neutrino 
interaction. At a large distance all neutrinos travel es- 
sentially coUinear so that |(v)r| —>■ 1, implying an ad- 
ditional "coUinearity suppression." A typical transverse 
component of the neutrino velocity is of order R/r if i? 
is the size of the source and r the distance. Therefore, 
a typical neutrino velocity along the radial direction is 
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FIG. 4: Radial evolution of the antineutrino polarization 
vector component for neutrinos emitted from a spheri- 
cal source. We compare the multi-angle solution (continuous 
curve that shows small oscillations) with the adiabatic limit 
of single-angle solutions where the dashed curve corresponds 
to Eq. (1331) and the dotted curve to Eq. pS)) . 



~ 1 — (i?/r)^, implying Cr cx r ^. Therefore, wc find 
the general scaling 



fir (X r 



(33) 



that is well known for spherical sources, and is also con- 
firmed numerically in Fig. [3) 

For a spherical source we have derived an explicit ex- 
pression for (vr) in Eq. ([M)) . In this case we find explicitly 



a = 4 



1- 



1 /i? 

2 v7 



(R/r) 



for 



- 1 



(34) 



Immediately at the source we find fiR = 3. The asymp- 
totic behavior for large r agrees with what one finds if 
one uses only a single angular bin with all neutrinos being 
radiated at 45° relative to the radial direction [2l|- 

In the papers by Duan et al., beginning with Ref. p^ . 
a somewhat different expression for Cr was used because 
they pictured the single- angle approximation as all modes 
oscillating like the radial one. Considering the average 
neutrino-neutrino interaction they find the equivalent of 



Cr = 



v^-{R/ry 



1 /R 
4 I r 



R/r 

2 



for 



(35) 



While the scaling at large radius is, of course, the same, 
their normalization is a factor of two smaller. 

We have checked with a numerical multi-angle exam- 
ple that indeed our result provides the proper approx- 
imation. In Fig. 2] we show the evolution of the an- 
tineutrino total polarization vector component from 



a multi-angle simulation (continuous curve). We com- 
pare it with the single-angle results corresponding to our 
approximation of Cr of Eq. (p4|) as a dashed curve and 
with Duan et al.'s prescription of Eq. ([55)1 as a dotted 
curve. In our numerical example, following Ref. [l^, 
we have fixed the neutrino-neutrino interaction strength 
at the neutrino sphere (radius R = 10 km) as /io = 
^/2GfFo = 0.7 X 10'"' km~^, the neutrino oscillation fre- 
quency ujQ ~ 0.3 km~^, the ratio of the number fluxes 
a ~ Fp^/F^^ = 0.75, and the neutrino mixing angle 
013 — 2.0 X 10~^. We have assumed the inverted mass hi- 
erarchy where bipolar conversions occur. For the single- 
angle cases, we have used the adiabatic approximation 
worked out in [3, H^j Ull , where depends on ra- 
dius only through /Jr and where therefore there are no 
nutation wiggles. We discuss in more details about the 
adiabatic approximation below in Sec. IIV CI 

Of course, the difference found between Duan et al.'s 
approximation and ours is not crucial in that the neu- 
trino fluxes and angular divergences assumed in previous 
studies had the character of toy examples. Since fi de- 
creases with r~'^, a factor 2 modification of /i translates 
into a 20% radial shift of fir that is not important for 
a toy model. The difference close to the source is much 
larger, but less important. Close to the source we typi- 
cally have strongly synchronized oscillations with a very 
small mixing angle, causing practically no difference in 
the overall solution. 

We finally consider the variation of fj. for a disk-like 
source. At a large distance, the variation with distance 
will be proportional to r""* as in the case of a spherical 
source where a factor of comes from the trivial ge- 
ometric flux dilution with distance, while another factor 
arises from the increasing collinearity of the neutrino 
trajectories with distance. To see this more explicitly we 
note that a typical "radial" velocity at a large distance is 
Vr = {1 — Vi — w|)^/^ where vi and V2 are typical trans- 
verse velocities in two orthogonal transverse directions. 
At a large distance we have wi <C 1 and U2 <C 1 so that 
the factor (l/|(v)xp - 1) m Eq. (I3l|) is w {vf + w|). A 
typical transverse velocity component of a neutrino at a 
distance r is R/r if R is the geometric size of the source, 
confirming the r"^ scaling of Cr at a large distance. 

Moreover, we can also derive the variation of /ix with 
zenith angle 8 at a large distance. If the source is a 
circular disk of radius R, the azimuthal component of 
the velocity, vi, will typically be vi cx {R/r). The 
polar component V2^ on the other hand, will vary as 
V2 oc (i?/r) cos 6, since the projected area of the source 
in this direction is suppressed by cos 6. This leads to 
Cr oc (1 + cos^ Q){R/rY. Since the flux itself varies as 
(i?/r)^ cos 8 we expect 



fir r 



cos 8(1 



cos^ 8) . 



(36) 



This is numerically confirmed and agrees with the varia- 
tion shown in Fig. [31 



10 



C. Adiabatic approximation 



V. COALESCING NEUTRON STARS 



The adiabatic approximation for collective neutrino os- 
cillations in the single- ang le case has been developed in 
a series of papers [3 [13 ■ If the neutrino flavor 

evolution is adiabatic, the flavor composition of neutri- 
nos as well as antineutrinos depends only on the initial 
state and the value of Wx//^x at each x, and there is no 
need to compute the evolution explicitly. 

Here, we follow and consider in particular the case 
of a system constituted by only two polarization vectors, 

P for neutrinos and P for antinuctrinos. At extremely 
large matter densities, the effective mixing angles are 

small, as a result P and P are aligned with each other 
and with B. In the limit of a large but slow varying neu- 
trino interaction strength fi^, the two polarization vec- 
tors move in a pure precession mode around B. During 
the evolution, the conservation of "flavor lepton number" 
implies 

Pz-Pz = l-a (37) 

is a conserved quantity, where a is the ratio of antincu- 
trino and neutrino fluxes. 

An explicit solution for the flavor evolution can be ob- 
tained simply in terms of a as follows. The results in 
(20| . when calculated in the limit of vanishing effective 
mixing angle, give 



1 — a + 2arj a (2 — a + ar\) 

2 VTT^ 



(1 — a + oa-])yJa(\ + r]) 
2y/2~~aTafj 



(38) 



where rj is the cosine of the angle between P and B at 
any location. It equals -1-1 to begin with, and varies be- 
tween -|-1 and —1 during the evolution. Eq. (j38p can be 
inverted to yield 77 at any location as a function of ti^x//^x- 

Knowing 77, one may compute the directions of P and P 
through 

Pz = ar] , 

P, = + (l-a) . (39) 

The values of P^ and Pz directly yield the survival prob- 
abilities Pee and Pee of Ve and Ve respectively. Pee and 
Pee are thus direct functions of the local value of lu^/ fix- 
In Fig. [4] we have shown how the analytical adiabatic 
approximation agrees very well with the multi-angle nu- 
merical simulation for the spherically symmetric SN case. 
In general, this result means that in the cases where neu- 
trino ensembles exhibit a self-maintained coherent be- 
haviour, the adiabatic approximation together with the 
single-angle limit described in Eq. (|3ip are sufhcicnt to 
describe the collective flavor evolution. 

On the other extreme, in the case of ensembles show- 
ing kinematical decoherence, the adiabatic approxima- 
tion is not appropriate, and the final outcome is a com- 
plete equalization of the neutrino fluxes. 



For purposes of illustration we now consider an explicit 
example for a disk-like source. Following the numerical 
studies of Ruffert and Janka [s^l for coalescing neutron 
stars, we use the following parameters for a disk source as 
an extremely schematic model. The all-species neutrino 
luminosity of the models shown by Ruffert and Janka 
varies between (1 — 12) x 10^^ erg s"^ and the luminos- 
ity in species other than i/e and De is always very small. 
Typically {E^^) < a typical average being 12 MeV 

that we will use for both species. In contrast to the SN 
case, here one expects an excess of the De number flux 
over Ve, so that in this case we normalize the polariza- 
tion vectors to the antineutrino flux. The ratio of number 
fluxes a = Fy^/Ffj^ varies between about 0.65 and 0.81 if 
we ignore one case where the number fluxes are exactly 
equal and one where the Ve flux practically vanishes. 

Inspired by these numbers and the pictures shown 
by Ruffert and Janka we adopt the following parame- 
ters for our homogeneous disk-like neutrino source: Ra- 
dius P 50 km, Lp^ = 3 X lO''^^ ^g s"!, {E^) = 12 McV, 
and a ~ 0.75. This implies a i?e flux at the disk of Pq ~ 
2.0 X lO''^^ cm~^ s~^ and a scale for the neutrino-neutrino 
interaction strength of //q = VSGp-Fo = 8.4 x 10~^ eV = 
4.3 X 10^ km~^. We have in mind two-flavor oscillations 
driven by the atmospheric neutrino mass difference for 



which we use Atti^^j^j 



2.4 X lO^'^ eV^ 



Then the aver- 
2 



age vacuum oscillation frequency is wq = = 
1.50 X 10"^° eV = 0.76 km~^. Here we have used 
that for an average neutrino energy (P) ~ 12 MeV 
and assuming a Maxwell-Boltzmann spectrum we have 
(1/P) = 1/(8 MeV). 

The transition between the synchronized and the bipo- 
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FIG. 5: Streamlines and surfaces of constant LJx//ix for neu- 
trinos emitted by an accretion disk with radius i? = 50 km. 
The surfaces are marked with the value of the survival 
probability, assuming the adiabatic limit of collective oscilla- 
tions. 
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lar oscillation region, where the flavor pendulum begins 
to nutate and neutrino transformations begin in earnest, 
occurs at [iBl 

i^.ii^, (40) 

which is 9 X 10"'^ for a = 0.75 in this scenario. 

In Fig. [5] we show the contours of constant ti^x/z^x, 
which are the same as those of constant Pee under adia- 
batic approximation, for inverted hierarchy. The contour 
with Pee = 1.0 is the one where the condition of Eq. (|40|) 
occurs and the synchronized to bipolar transition takes 
place. Within the region enclosed by this contour, the os- 
cillations are strongly synchronized so that macroscopi- 
cally nothing much happens. This contour delineates the 
area where the flavor survival probability is essentially 
unity for both Vf. and Vf^. Outside this region, neutrino 
transitions begin and we show contours of the Ve survival 
probability that indicates to which extent the initial 
and Vf, fluxes have been converted to other flavors. 

For our specific numerical example, we find that col- 
lective neutrino flavor conversions start at z > 400 km. 
Since the v^Ve annihilation rate per unit volume decreases 
very rapidly far from the source (for a spherical geometry 
one would expect ~ [i^l) we do not expect that these 
flavor transitions affect significantly the neutrino energy 
deposition rate in the plasma. However, the features of 
neutrino fluxes emitted from coalescing neutron stars are 
rather model dependent. Here we speculate that if the 
asymmetry between the emitted neutrino species would 
be smaller, one could have collective flavor conversion 
or perhaps some sort of kinematical decoherence of the 
neutrino ensemble in a region close to the source, with a 
possible impact on the annihilation rate. We also note 
that the matter density in the region close to the source 
is expected to be p ~ lO^""^'^ g cm^'^, thus it is so strong 
to prevent ordinary resonant flavor conversions [4l| . At 
larger distance one expects that two resonant level cross- 
ings will occur. However, due to the uncertainties of the 
matter density profile in these regions, it is difficult to 
predict their effects on the neutrino burst. 

VI. CONCLUSIONS 

Neutrinos streaming from powerful astrophysical sour- 
ces such as SN cores or coalescing neutron stars are so 
dense near the source that they must show nonlinear fla- 
vor oscillations induced by the neutrino-neutrino refrac- 
tive effect. Numerical simulations reveal a rich variety 
of phenomena, some of which have been explained with 
simple analytic models. However, numerical simulations 
thus far have been restricted to homogeneous gases evolv- 
ing in time or to sources with exact spherical symmetry. 
More general geometries are numerically much more de- 
manding and have not yet been studied. 

Therefore, we have studied what might be expected un- 
der the assumption that the multi-angle instability plays 



no role and that the neutrino ensemble is largely char- 
acterized by self-maintained coherence. In this case one 
is led to a unique formulation of the collective equations 
of motion that imply that collective flavor oscillations 
should be thought of as a one-dimensional phenomenon 
along the streamlines of the underlying neutrino flux. 
Close to the source these streamlines are usually curved 
even though, of course, the underlying neutrino trajec- 
tories are straight. (We have neglected the gravitational 
bending of trajectories.) Therefore, even if the neutrino 
stream has no global symmetries, the collective oscilla- 
tion problem is relatively simple. 

We have used the concept of "self-maintained coher- 
ence" in the most restrictive sense that applies when 
the neutrino gas is dense, i.e., when a typical neutrino- 
neutrino interaction energy /i is large compared to a typ- 
ical vacuum oscillation frequency wo = {Am^ /2E). The 
neutrino ensemble in this case evolves along a stream- 
line as one unit that previously has been described as a 
gyroscopic pendulum in flavor space. All neutrino and 
anti-neutrino polarization vectors point essentially in the 
same direction in flavor space, the pendulum direction, 
allowing for the simplifications that lead to our collective 
equations. We have provided a prescription for defin- 
ing the effective neutrino-neutrino interaction strength fi 
that works for general source geometries. Our result is 
somewhat different from what was used in the previous 
literature, but our expression for fj, agrees better with 
numerical simulations. 

There are more general forms of collective behavior. In 
a homogeneous isotropic system with a density that de- 
creases adiabatically in time, the polarization vectors are 
at first aligned and stay in a single co-rotating plane even 
when n becomes of order wq and smaller. However, they 
will align or anti-align themselves in the mass direction, 
leading to the phenomenon of spectral splits. This is a 
more general case of self-maintained coherence beyond 
our simple picture and is analytically fully understood. 
In the spherically symmetric case, numerical multi-angle 
simulations show a similar behavior. If the asymmetry 
between neutrinos and antineutrinos is large enough to 
prevent multi-angle decoherence, the neutrino and an- 
tineutrino modes evolve collectively essentially along the 
co-rotating plane, but with intriguing three-dimensional 
patterns that evolve collectively and are stable in flavor 
space. Once more spectral splits develop, similar to the 
single-angle approximation. 

We speculate that these more general forms of behavior 
also occur for more general source geometries if we follow 
the neutrino stream lines. It is clear that multi-angle de- 
coherence will occur if the asymmetry between neutrinos 
and anti-neutrinos is extremely small. Even a few per 
cent asymmetry, on the other hand, implies that close to 
the source we have synchronized oscillations, or, in the 
pendulum picture, it precesses fast without nutations. 
Therefore, close to the source where the streamlines are 
curved, our treatment should be most appropriate. If the 
effective mixing angle is small (as it will be because of the 
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presence of ordinary matter) , superficially no flavor oscil- 
lations happen close to the source and the "true action" 
begins at the synchronization radius where the inverted 
pendulum begins to nutate and begins to lose its initial 
orientation. If this point is reached at some distance, the 
streamlines will essentially point straight away from the 
source and naively one should think that henceforth the 
evolution is similar to the spherically symmetric case. In 
fact in the adiabatic approximation, one only needs to 
know the local value of Wq/m in order to determine the 
neutrino survival probabilities, for any source geometry. 

Granting these assumptions, we have briefly studied a 
toy model for coalescing neutron stars where the source 
is taken to be a homogeneously radiating disk. For our 
numerical example, we find that the region where collec- 
tive oscillations will convert the original and fluxes 
is quite far from the emission surface. Therefore these 
flavor conversions should have a negligible effect on the 
neutrino energy deposition in the plasma. However, due 
to the uncertainties of the original neutrino fluxes, we can 
not exclude situations with a smaller asymmetry between 
the original fluxes. Such tiny neutrino asymmetries may 
give rise to multi-angle decoherence, that could produce 
signiflcant flavor conversions near the emitting disk. This 
question ultimately will need to be addressed numerically 
by means of large-scale numerical simulations of neutrino 
flavor evolution near a non-spherical source, in analogy 
to what has been done in the case of neutrinos streaming 
off a SN core. 

Our study suggests that self-maintained coherence 
among different neutrino modes may well be a typical 
form of behavior even in non-spherical systems and no- 
tably in the interesting case of coalescing neutron stars. 
The flnal verdict on the role of collective neutrino oscil- 
lations in such systems can only come from numerical 
studies. At very least, our approach provides a simple 
limiting case against which one can measure the output 
of numerical simulations. 



APPENDIX A: DISK SOURCE 

We here briefly describe the numerical integrations for 
a disk source that were used in the main text. The direc- 
tion of neutrino momenta are determined by two polar 
angles Op and (j)p, by which we can express the compo- 
nent of the neutrino velocity in x and z direction respec- 
tively by 



V 



cos 
cos 0< 



<p sin Op , 



p ■ 



(Al) 



Our flrst task is to determine the flux Fx, i.e., the x 
and z components F^ .^ and F^^^, each as a function of x 
and z. These are given by 



1 

8^ 
1 

8^ 



dE d cos t 



= zz;^ I ^Ed cos Op, 



'!>pw^/x,z(w'',w^), 



One can also define the local neutrino density 

nx,z = -^ J dEdcosOpd(j)pf^Av''-,v'') . (A3) 

In the case of synchronization, from Eq. (j2ip we arrive 
at the following EOM 



(A4) 



where 



X.Z 



F^ 



(A5) 



In this case, we have obtained a first order partial dif- 
ferential equation (PDF). Since {v^ ^ and ^ are not 
constant, in this case the streamlines will be curved. If 
the streamlines are graphs of the hmction z{x), it follows 



dz 
dx 



(A6) 
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(i.e., the tangent line to the graph of z{x) at (x,z) is 
parallel to v.) The ordinary differential equation (ODE) 
in Eq. (|A6p is the so-called characteristic equation for the 
associated PDF. Its solutions are the streamlines for the 
PDF. 

The range of the previous angular integrations are de- 
termined by the coordinates x and z. It could be conve- 
nient to express the previous integrals at a given point, 
writing Op and in terms of the coordinates tq and 
00 on the disk, by tracing back the direction of neu- 
trino momenta on the disk, as shown in Fig. [B] Let 
the neutrino emitted from a point (ro cos <^oi ''o sin^o, 0) 
pass through the location (x,0, z). The direction of 
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so that 





FIG. 6: Relation between the momenta variables 9p and (j>j, 
and coordinates ro and 4>o on the disk. 

the momentum of this neutrino is along the unit vector 
(sin 9p cos 4>p, sin 9p sin 0p, cos Op). Therefore, 

{x — ro cos(/)o, — 7'o sin(/)o, z) cx 

(sin 0p cos (/)p, sin^p sin 0p, cos 0p) . (A7) 

Wc can now write down Op and (j}p explicitly in terms of 

tan Op = ^ri+x--2xr,cos^ 
z 



SO that 



COS^?p 



sin^p 



For the angle t^p, we get 



i/rp + — 2xrQ cos 0o + z'^ 

^Jrl+ x'^ - 2xro cos^ ^^^^ 
-^rp + — 2xro cos 0o + 



sin0p 
cos 4>p 



rp sin (/io _ rp sin (/>o 



z tan 6'p 
ro cos ( 



sinflp 

a- - ro cos 4)p _ x ~ rp cos (AlO) 
sin Op z tan 0p ' 



4>p = ~ arctan 



ro sin 0o 



(All) 



^ X — rp cos (pp ^ 
Using as variables of integration ro and ^o, we obtain 
dE 



dcosOpd(j)pf^,^{v'',v'') 



drod</)odetJ(ro,</>o)/.,z(«",t'") , (A12) 



where J{rp,<f>p) is the Jacobian that relates the two co- 
ordinate systems, and < ro < 1 and < (f>o < 27t. 

The derivatives relevant for the Jacobian are the fol- 
lowing 



d cos Or, 



-z{rp — X cos 4>p) 



drp (rp + a;2 _ 2j;ro cos ^p + 2;2)3/2 



d cos Op 
d(j)o 



'ZxrQ sin 0o 



Uipp 

drp 
d(j)p 



(rp + x"^ ~ 2xrp cos 4)p + z2-)3/2 
a; sin 0o 
rg -I- a;2 _ 2a;ro cos (jip ' 
rQ — xro cos 4>p 



hp r^+x^ — 2xrp cos (j>p 



(A13) 



so that 



^ dcosOpd4>p d cos Op d4>p 

drp d(j)p d(j)p drp 

_ ~zrp 

(rg + x'^ — 2xrp cos <j)p + z^)'^!'^ 



(AM) 



The net result is 



Ft 



<^E -zrp{x-rpCOS(t)p)fx,z{v'^-,v'') 
— d^odro ^ _ _ ^ ^^^^ 



87r3 
dE 
8^' 



^;--^d0odro7;3 



(rp -I- _ 2a:;ro cos (fip + z'^Y 
z{-zrp)fx^;,{v'',v'') 



(ro + x'^ — 2xrp cos (?f)o -I- 2:^)^ 

(A15) 
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